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Abstract. - The functional RG (FRG) approach to pinning of d- dimensional manifolds is 
reexamined at any temperature T. A simple relation between the coupling function R(u) and 
a physical observable is shown in any d. In d — its beta function is displayed to a high order, 
ambiguities resolved; for random field disorder (Sinai model) we obtain exactly the T — 
fixed point R(u) as well as its thermal boundary layer (TBL) form (i.e. for u ~ T) at T > 0. 
Connection between FRG in d — and decaying Burgers is discussed. An exact solution to the 
functional RG hierarchy in the TBL is obtained for any d and related to droplet probabilities. 



Elastic manifolds pinned by quenched disorder [1] are the simplest system to study glass 
phases where (dimensionless) temperature is formally irrelevant, scaling as TV = TL~ with 
system size L. They are parameterized by a displacement (Af-component height) field u(x) = 
u x , where x spans a d-dimensional internal space. The competition between elasticity and 
disorder produces rough ground states with sample averages (u x — u x ') 2 ~ \x — x'\ 2( > (6 = 
d — 2 + 2£). These are believed to be statistically scale invariant, hence should be described 
by a critical (continuum) field theory (FT) . The latter seems highly unconventional in several 
respects. First, an infinite number of operators become marginal simultaneously in d = 4 — e. 
This is handled via Functional RG methods where the relevant coupling constant becomes a 
function of the field, R(u), interpreted as the (second cumulant) disorder correlator [2]. A more 
formidable difficulty then arises: at T — both R(u) and, more generally, the full effective 
action functional T[u], appear to be non-analytic ( 1 ) around u — 0. A linear cusp in R" (u) was 
found in one loop and large A^ calculations [2,3]. Qualitative (two mode minimization) and 
mean field arguments relate this cusp to multiple metastable states and shock type singularities 
in the energy landscape [4]. As a consequence, ambiguities arise in loop corrections [5]. 
Although candidate renormalizable FTs have been identified [5,6] (working directly at T = 0) 
this problem has, until now, hampered derivation of the field theory from first principles (with 
the notable exception of the N = 1 depinning transition [7]). 

Working at non-zero temperature T > should help define the theory, and T[u] has been 
argued to remain smooth within a "thermal boundary layer" (TBL) of width u ~ Tl around 
u = 0. This width however shrinks as Tl — > in the thermodynamic limit, and if a fixed u, 
large L limit exists for any fixed small T it should unambiguously define the (non-analytic) 

( 1 )the non-analyticity of T occurs at finite scale (the Larkin length) contrarily to e.g. the critical 4 theory 
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"zero temperature theory" : this program, called " matching" , was proposed and extensively 
studied in Ref. [8]. It does, to some extent, rely on a scaling ansatz proposed there for the 
TBL. This ansatz was shown to be consistent to 1-loop with the droplet picture [8] and, in the 
(near) equilibrium (driven) dynamics, to account for the phenomenology of ultraslow activated 
(creep) motion [8,9]. Although its physics is reasonable, it is not yet established how a critical 
renormalizable FT emerges from it as Tl — > 0, with a finite unambiguous beta function. 

Another field of physics where an (unconventional) field theoretic description is needed, but 
remains elusive, is high Reynolds number turbulence. There too the scale invariant regime, 
the inertial range, needs regularization at the small dissipation scale set by the (formally 
irrelevant) viscosity v [10]. Connections between these two tantalizing problems can be made 
quantitative within the simplified Burgers turbulence, a much studied problem [11-14]. 

Given its central role in the FRG, it is of high interest to obtain the precise physical 
content of the (fixed point) function i?(u), beyond previous qualitative arguments. In the 
FT, a precise, but abstract, definition was given, from the replicated effective action at zero 
momentum, which allowed for a systematic dimensional expansion. From it, it was observed 
that R"(0) gives the exact sample to sample variance of the center of mass of the manifold 
(a typical observable with a universal T — limit), while R""(0) yields sample to sample 
susceptibility fluctuations (a finite temperature observable which diverges as Tl — ► 0). It 
would be useful to relate directly the full function R(u) to an observable and cleanly separate 
zero from finite T contributions. 

The present Letter is a short account of a recent study [15] aimed at clarifying the physics 
encoded in the FRG and its connections to Burgers turbulence. We obtain a simple operational 
definition valid in any d, not only for R(u), but also for higher cumulants, and the full 
(replicated) effective action T[u]. It makes explicit its T = physics and at T > makes 
precise the relation between the TBL form of the effective action, droplet probabilities, and 
dilute (functional) shocks, via a (functional) decaying Burgers equation. Next, the instructive 
d = case is studied. For N = 1, the matching program started in Ref. [8] is pushed 
to obtain here the (unambiguous) beta function to four loop, and related to works on the 
inviscid distributional limit of Burgers equation [11]. In the sub-case of the Sinai (i.e. random 
field) model, the exact R(u) is computed at T = 0. The TBL rounding form at T > is 
also obtained. Obtaining the thermal rounding form in any d amounts to solve an infinite 
hierarchy of (functional) exact RG equations: remarkably, this can be achieved, the solution 
being parameterized by droplet probability data. All details are given in [15]. 

The model studied here is defined by the total energy: 

H V[U] = ~ j 9xy u xUy+ j V(U X ,X) (1) 
^ J xy J x 

in a given sample (u G R N ). The distribution of the random potential is translationally 
invariant, with second cumulant V(u,x)V(u',x') — S d (x — x')Rq(u — u') and V(u,x) = 0. 
This implies the statistical (tilt) symmetry (STS) under (x,u x ) — ► (x,u x + 4> x ) [1]. Several 
results here are valid for arbitrary g xy , but we often specialize to g^ 1 = q 2 + m 2 in Fourier 
space, where the small mass provides a confining parabolic potential and a convenient infrared 
cutoff at large scale L m = 1/m. In all formula below one can replace f = j d d x — ► J2 X 
and 6(x — x') — > 5 xx i in the bare disorder correlator, i.e. a lattice provides a UV cutoff 
which preserves STS Numerous physical systems are modelled by e.g.: (i) CDW or vortex 
lattices [1], for a periodic Rq{u), (ii) magnetic interfaces with bond disorder, for a short range 
(SR) Rq(u), (iii) magnetic interfaces with random field disorder, of variance a, for a long 
range Rq(u) ~ ~cr|u| at large u. 
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Let us briefly recall the convenient definition of R(u) used in the FT. The model is studied 
using replica fields u%, a = 1, ..,p, with bare action: 

S[u] = ^ E / 9*v« ~ VFi E / R o« <) ( 2 ) 

ZJ a Jx V L± ab Jx 

and disorder-averaged correlations of (JTJ identify with replica correlation functions of at 
p = 0. These (the connected ones) are obtained from Taylor expanding the W functional, 

From it one defines, via a Legendre Transform, the 
effective action of the replica theory, T[u] = J x J2 a u xjx ~ W\j]. It generates (in a Taylor 
expansion in u) the renormalized vertices, i.e. those where loops have been integrated, and is 
thus the important functional for the FRG. To define the renormalized disorder one assumes 
an expansion in number of replica sums: 

rr«i - V [ 9 *y u> y V R[uab] V ^-1 + ™ 

a xy ab abc 

where STS implies that the single-replica term is the bare one, and the form of the n-replica 
terms, e.g. R[u ab ) is a functional depending only on the field = u% — u h x , whose value for 
a uniform field (i.e. local part) defines R(u), i.e = u}} = L d R(u) ( 2 ). It was used in 

the FT [5,6] to compute the beta function, — md m \n R(u) = /3[R](u), in powers of R, and its 
derivatives. 

We have shown that this abstract definition is equivalent to a physical one: for each 
realization of the random potential V, one defines the renormalized potential functional V[v] = 
t^[{t>a}] as the free energy of the system when centering the quadratic potential around u x = 

V[{v x }} = -Tin f Y[du x cM-HvM/ T ) ( 4 ) 

*' X 

Hy.vH = \ I 9xy( u x ~ V x )(u y ~Vy)+ { V{U X ,X) 
J xy J x 

Using STS one sees that the renormalized energy landscape has second cumulant correlations: 



V[{vx}]V[{v' x }}=R[{v x -v' x }} (5) 

and V = (averages are w.r.t. V). The result shown in [15] is that R = R. Hence one can 
measure the 2-replica part of the effective action by computing the free energy in a well whose 
position is varied. Choosing a uniform v x = v, one obtains its local part: 



V(v)V(v') = L d R(v - v 1 ) (6) 

where V(v) — V[{v x — v}], using a parabolic potential centered at u x — v. Performing 
the Legendre transform [15] (more involved) relations are found for higher cumulants, e.g. 
S = S - 3symabcgxyR' x [v ab }R' y [v ac }. The STS property was used: for a non STS model, 
e.g. with discrete u, either it flows to the STS fixed point as m — > 0, and the above holds 
asymptotically, or it does not and a (more involved) extension holds [15]. 

( 2 )such zero momentum renormalization conditions are standard in a massive theory. It is not presently 
known how to close FRG using other conditions, e.g. symmetric external momenta, as in massless theories 
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From Q the renormalized (pinning) force functional T x = V£[v] = SV[v]/Sv x is related to 
the thermally averaged position in presence of the shifted well, via T x = J g~y (v y — (u v )h v „ ) • 
Hence the force correlator functional R xy [v] has a nice expression. For uniform v x — v and 
at T — it is simple: denote u x (v) the minimum energy configuration of £Zy>M for a fixed 
v x = v and u(v) — L~ d f u x (v) its center of mass position. Then, denoting A(v) = —R"(v): 

(v - u(v))(v' - u(v')) = A(v - v')L- d m- i (7) 

which generalizes to non-zero T (replacing u(v) by its thermal average) and to the full multi- 
local functional ( 3 ) 

For fixed L/a, where a is the UV cutoff scale, the minimum is expected unique for con- 
tinuous distributions of V , except for a discrete set of values v„ which are positions of shocks 
where u(y) switches between different values (e.g. u\ to U2) and the force is discontinuous (at 
T = 0): below, the strength of each shock is noted ttg* = Ui — u\. 

The renormalized pinning force satisfies an exact RG (ERG) equation (with dg = —md m g): 

- 2md m T x M - f dg yz (T^ yz [v] - T' xy [v]T z [v] ) (8) 

J yz 

a functional generalization of the decaying Burgers equation to which it reduces for d = 0: 

d t F(v) = |f» - F'(v)F(v) (9) 

with t — m~ 2 , F(v) — V'(v), usually written d t u + u' x u = ^u xx , identifying u, x, v in Burgers to 
F,v,T/2 in the FRG (while V[v] satisfies a functional KPZ-type equation). The stochasticity 
in (JHJl,© comes from their (random) initial conditions F(v) = V'(v) and T x [v] = V x [v] for 
t = 0, m = 00. Eq. Q (and its primitive) is equivalent to an infinite ERG hierarchy for the 

n-th moments S , (")(«i,2,..,n) = (-) n V r (vi)..T>(w„) in d = 0: 

2T 2 - 

-md m R(v) = -^R"{v) + -^5uo(0, 0, v) (10) 

-«4S W («l^) = ^[S££ («l,2,..»)] + ^[ffi.0 ) («l,l,2..n)] (11) 

where S = = S, subscripts denote partial derivatives and [..] is symmetrization. A 

similar, more formidable looking functional hierarchy exists for any d: 

- md m R[v] = Tdg xy R%[v] + dg zz ,SH?[0, 0, v] (12) 

together with ERG equations for S and higher moments. In both cases a related hierarchy 
exists for the cumulants R, S, .. defining T[u] in studied in [8,16]. The usual RG strategy 
is to truncate them to a given order in R yielding the beta function. Ambiguities in the limit 
of coinciding arguments in 111I12|I may arise in doing so directly at T = 0. 

We start with d = (and N = 1), a particle in a ID random potential V(u), aiming 
to obtain an unambiguous beta function as m — > 0. We define rescaled T = 2Tm e and 
R(u) = im e_4 ^i?(um^) [this should yield a FP when correlations of V grow as u 9 It]. Trying 
first standard loop expansion at T > (R analytic), we obtained from (jl 1|> the beta function 

Q generalize to any N, and to two copies as used in chaos studies [17] Vi(v)Vj(v') = L d Rij{v — v'). 
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Fig. 1 - Construction of the joint probability P({Ei,Vi}) that V(vi) = Ei at points v%: the random 
walk V(u) must remain above all parabola centered on the Vi of apex Ei intersecting at points u* . 
Each independent interval [u*,u* +1 ] can be treated as in [18]. 



— md m R\n = (3[R,T]. To n loop, it is a sum of terms of order T p R n+1 p , < p < n. 
The one-loop equation (i.e. adding TR" to the first three terms in (|13|l below) exhibits the 
standard TBL for u ~ f discussed in [8]. To 2-loop a term ^TR""(Q)R"(u) appears, and 
using the TBL identity lim m _>o TR""(0) = i?"'(0 + ) 2 , exact at one loop, produces precisely the 
2-loop "anomalous" term in l|13f) below. Alas, one finds [15] that this procedure fails at 3-loop. 
One must instead examine the whole ERG hierarchy as in [8] . There, a method to obtain the 
unambiguous beta function was found by verifying order by order, a continuity property of the 
T-cumulants S^" i(ui.. n ) upon bringing points together. We completed in [15] the derivation 
of the (local) beta function, obtaining (up to a constant, with R" = R" — R"(0)): 

-ma m R - (« - 40* + Cug + - «"(0)fi»] + \ {{ R"r- - «"'(0+) 2 )R" (13) 

+ Lm"f(i("'f + —((R'"f - R"'(0+ ff + -R"((R"'fR"" - B"'(0+) 2 B""(0 + )) 
16 32 4 

+ J_(R")3(^(5))2 + ^( R "fR"'R""R^ + -R"((R"') 3 R^ - R"'(0+) 3 R ( - 5 h0+)) 
96 16 8 

+ ^(R")2(^"")3 + i R «((fl'")2(fl"'')2 _ li?"'(0+) 2 (i?"") 2 - -R"'<0+) 2 R""<0+) 2 ) 

16 16 6 6 

+^((^"') 2 ~ R"'(0 + ) 2 ){{R"') 2 R"" + -^R""R"'{0+) 2 - ^i?""(0+) J R"'(0+) 2 ) + 0(R 6 ) 

The first line are one and 2-loop terms, the second is 3-loop, the last three are 4-loop. Normal 
terms (i.e. non vanishing for analytic R(u)) are grouped with anomalous "counterparts" 
to show the absence of 0(u) term, a strong constraint (linear cusp, no supercusp): these 
combinations can hardly be guessed beyond 3 loop. This shows the difficulty in constructing 
the FT, already in d = 0. We emphasize that (|13(l results from a first principle derivation. 

R(u) being a physical observable, we look for cases where it can be computed. The 
Brownian landscape V(u), the so-called Sinai model, is interesting as the d = limit of random 
field disorder. Recently we obtained the full statistics of (deep) extrema in presence of a 
harmonic well [18]. This is generalized [15] as described in Fig^ Graphically the renormalized 
landscape V(v) = E is constructed by raising a parabola P v , y(u) = — ^-(m — v) 2 + E' from 
E' = — oo until it touches (for E' = E) the curve y = V(u) at point u = m(v), position 
of the minimum of Hy iV (u), E being the maximum (apex) of the parabola. P v touching at 
two points ui(v) < v,2(v) signals a shock at u = v. Computing P(E\, V\; E 2 , V2) (see Fig. 
yields [15] (for m 2 = 1, a = 1) and R(v) = R(v) - R{0): 
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!x 1 Jx 2 Ai{i\ 1 )Ai(i\ 2 ) Ai{iX x )Ai{iX 2 ) 

(14) 

where a — 2 -1 / 3 , ba 2 = 1, J x = ^ and all integrals converge well. One finds R(v) « 

-v + 0.810775 at large v, and recovers [18] v? = -i?"(0) = 1.054238. Once rescaled, lfH)l 
should be a FP of (|T3|) corresponding to £ = Crf =4/3. 

At non-zero T, one reexamines (0} taking into account, within a droplet calculation, the 
probability density D(y)dy for two degenerate minima of V, spatially separated by y — u 2 — u\ 
(D(y) = D{~y)). It yields the TBL form (for v ~ Tm 2 ): 

R"{v) = R"(Q) + m 4 T(y 2 F 2 (m 2 yv/T)) y (15) 

with F 2 (z) = | coth | — i and (..),, = J dy..D(y) is normalized by the STS identity (y 2 ) y = 
2 /in 2 . Since F 2 (z) ~ |z|/4 at large z (f 1 51) yields consistent matching between finite T (droplet) 

quantities in the TBL and the cusp of the T = FP for v = 0(1), with R'"(0+) = ^^f- 
(|15|l should be more generally valid in d = (any N), but in the RF case it is known [18] 
that (setting m = a = 1) D(y) = i J x ^ x ^ Al Ai(i\ 1 )Ai 2 (i\2) — found to be consistent with 
i?"'(0+) = 0.901289 from 12} . Remarkably, generalizes to higher moments S^ n \ with 
generalized functions F n explicitly obtained in [15] yielding an exact "droplet" solution of the 
hierarchy {TTJ within the TBL. 

Since in d = the FRG I© identifies with decaying Burgers, we emphasize the correspon- 
dence: 

- R"(0) = u(zp , ^i?""(0) = v{Vu{x)) 2 = I (16) 



(more generally u(x)u(0) = — i?"(x)), both have finite limits as v — > 0. The second is the 
dissipative anomaly: also present in 3D Navier Stokes. In Burgers it is due to shocks. The 
(equivalent) finite limit of the l.h.s. implies a thermal boundary layer in the FRG. Dilute 
shocks in Burgers are equivalent to droplets and a TBL in the FRG where u 2 i == u(0 + ) — 
u(0~), and ifFST) can indeed be recovered from a single shock solution in Burgers upon averaging 
over its position [15]. The celebrated Kolmogorov law in the inertial range: 

1 



2 S m (0,0, M ) ~ eu = — (u(x) - u(0))3 ~ -ex (17) 

corresponds to the non-analytic behaviour of the third cumulant at small argument in the 
T = theory. Identical coefficients in (|16|l and (|17J) are a consequence of matching across the 
TBL (i.e. viscous layer), identifying the second derivative of l|10(l at v = (for v > 0) and 
v = 0+ (for v -> 0) i.e d t R"(0) = TR""{Q) = d t R"(0 + ) = S 112 (0, 0, 0+). Similar relations 
exist in stirred Burgers (and Navier Stokes) [10]: there the dissipation rate 1 is balanced by 
forcing, instead of scale-invariant time decay of correlations, but small-scale shock properties 
should be rather similar. Closure of hierarchies similar to (jl 1|> was proposed there [13] in terms 
of an "operator product expansion". Recent studies cast doubt on such simple closures [11]: 
N=l decaying Burgers (and stirred [14]) can be constructed in the inviscid limit {y — > 0) 

using distributions, e.g. tF'(v) = 1 — J2 s u< 2i ~ v s)- It is shown there that shock "form 
factors" (i.e. size-distribution) determines small distance (non-analytic) behaviour of moments 

of velocity differences, (F(v) — F(0)) p ~ [i p vsign(v) p+1 , with [i p = X^ s ( u 2i ) P H V ~ v s)- In the 
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FRG these are equivalent to droplet distributions as we show [15] that /i p = (\y\ p+1 ) y /(y 2 ) y , 
e.g. consistent with R'"(0+) = u(0+)Vu(0) = fi2/(2t 2 ) given above. The T = distributional 
limit of (jHl derived in [11] is equivalent to dtV(v) — —^F(v + )F(v): it validates the first- 
principle FRG discussed above yielding l|13|) (the central property being continuity of all S1..1 
since F(v) remains bounded). These considerations should be universal for dilute shocks, i.e. 
independent of details of shock probability distributions. For RF disorder the full distribution 
of shock parameters {u 2 i,v s } is known exactly [19]. It is used in [15] to obtain from Q) 
another expression for A(t>) fully consistent with (|14|) . 

Extensions to higher d are studied in [15]. Let us indicate here that similar droplet esti- 
mates can be performed and yield an exact solution of the full functional hierarchy (|12|) for 
all moments. Within the TBL m 2 v x /T — 0(1) the R[v] functional reads: 



where H'^z) = F2(z). To find this solution one considers a small density (of order Tm e ) 
of well separated "elementary droplets", i.e. local GS degeneracies u(l x = — u\ x . (..)d 
denotes the average over them. Eq. I|18f) relates droplet probabilities to the TBL in the FRG. 

To conclude we related FRG functions, e.g. R(u), to observables. This allows to compute 
them in simple cases, and provides a method to measure them in numerics and experiments. 
Their relations to shocks in energy (or force) landscape was made precise, via a generalized 
Burgers equation. We have shown how shock form factors and droplet distributions are 
related to the FRG functions. Questions such as the extent of universality in the TBL, how 
do properties of N = 1, d = Burgers extend to functional shocks (e.g. Kolmogorov law) 
remain tantalizing but can now be adressed. 

We warmly thank K. Wiese and L. Balcnts for enlightening discussions and long standing 
collaborations on FRG. We acknowledge support from ANR under program 05-BLAN-0099-01 
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